Athermal energy loss from X-rays deposited in thin superconducting films on solid 

substrates 


Alexander G. Kozorezov 1 , Colin J. Lambert 1 , Simon R. Bandler 2 , Manuel A. Balvin 2 , Sarah E. Busch 2 , 
Peter N. Nagler 2,3 , Jan-Patrick Porst 2,3 , Stephen J. Smith 2 , Thomas R. Stevenson 2 , John E. Sadleir 2 
1 Department of Physics, Lancaster University, Lancaster, UK, 

2 NASA/Goddard Space Flight Center, Greenbelt, MD 20771, USA, 

3 Physics Department, Brown University, Providence, RI 02912, USA 
(Dated: February 5, 2013) 

When energy is deposited in a thin-film cryogenic detector, such as from the absorption of an X-ray, 
an important feature that determines the energy resolution is the amount of athermal energy that can 
be lost to the heat bath prior to the elementary excitation systems coming into thermal equilibrium. 

This form of energy loss will be position-dependent and therefore can limit the detector energy 
resolution. An understanding of the physical processes that occur when elementary excitations are 
generated in metal films on dielectric substrates is important for the design and optimization of a 
number of different types of low temperature detector. We have measured the total energy loss in 
one relatively simple geometry that allows us to study these processes and compare measurements 
with calculation based upon a model for the various different processes. We have modeled the 
athermal phonon energy loss in this device by finding an evolving phonon distribution function that 
solves the system of kinetic equations for the interacting system of electrons and phonons. Using 
measurements of device parameters such as the Debye energy and the thermal diffusivity we have 
calculated the expected energy loss from this detector geometry, and also the position-dependent 
variation of this loss. We have also calculated the predicted impact on measured spectral line-shapes, 
and shown that they agree well with measurements. In addition, we have tested this model by using 
it to predict the performance of a number of other types of detector with different geometries, where 
good agreement is also found. 


I. INTRODUCTION 

An understanding of the key processes that take 
place when energy is deposited is important for any type 
of radiation detector. In this paper we present direct 
measurements and detailed calculations of energy loss 
and fluctuations in this loss from athermal phonons that 
escape from a metal film into a dielectric substrate. This 
loss mechanism can influence the measured energy reso- 
lution performance of a wide variety of particle detectors. 

There are two major classes of low temperatures 
photon detectors. The first is a class of devices operating 
under non-equilibrium conditions; the best examples of 
these are superconducting tunnel junctions (STJs) and 
microwave kinetic inductance detectors (MKIDs). For 
these detectors non-equilibrium quasiparticles are cre- 
ated in the process of energy down-conversion which re- 
sult in the breaking of numerous Cooper pairs. This 
down-conversion process typically lasts from several tens 
to hundreds of picoseconds. During the down-conversion 
process the photon energy is transformed into elec- 
tronic excitations and phonons, the properties of non- 
equilibrium distribution of quasiparticles and phonons 
being used to determine the energy deposited. The mech- 
anisms for determining the energy partition between the 
two systems far away from equilibrium are important and 
lead to the well-known statistical Fano-limit for the res- 
olution of this type of detector. 

A second class of low-temperature detectors are 
equilibrium devices. These include microcalorimeters 
and bolometers that are used to measure the energy of 


single photons and the power from a flux of photons re- 
spectively. In these devices the different types of exci- 
tations rapidly come into thermal quasi-equilibrium at a 
temperature slightly above that of the thermal bath. We 
consider the distributions to have reached quasi-thermal 
equilibrium when the electron and phonon systems are in 
equilibrium with each other, but the temperature of these 
systems is still changing. The measured signal is deter- 
mined from the mean energy of all excitations and is de- 
termined by the temperature excursion. In low tempera- 
ture (<1K) microcalorimeters thermal quasi-equilibrium 
is typically established on the time-scale of less than 1 fj, s, 
and Fano- fluctuations are completely washed-out. 

One key question that is relevant to all types of de- 
tectors is whether the whole deposited energy is retained 
in the sensing element and converted into the response 
signal. Energy that is leaked from the sensor prior to 
generating the ultimate number of quasiparticles or es- 
tablishing thermal quasi-equilibrium may naturally oc- 
cur in all types of detectors. Although the energy down- 
conversion process is usually fast, a significant fraction 
of the photon energy can be accumulated within the en- 
ergetic phonon system (athermal phonons). Depending 
upon the geometry of the detector, these phonons have 
some probability of escape from metal films into a di- 
electric support, such as a silicon wafer. This can result 
in a significant deficit in the energy detected. A num- 
ber of different detector types have athermal phonons 
produced in films fabricated on solid substrates, where 
athermal phonon escape can degrade performance either 
from the deficit in the amount of energy finally detected, 
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or through variations in the amount of energy lost de- 
pending upon the location of the initial photon interac- 
tion. An example of performance being affected by the 
deficit in energy is in single photon detectors operating 
at visible and near-infrared wavelengths 1 . Here the spec- 
tral resolving power is typically not very high (~10), and 
therefore variations in energy loss ~1% are not impor- 
tant. But the loss in energy, as high as 50% in these low 
energy photon detectors, can significantly affect the per- 
formance. An example where the broadening of energy 
resolution is more important than the energy deficit is in 
some X-ray detectors fabricated on solid substrates with 
STJ sensors. Variation in athermal phonon escape was 
identified as the mechanism responsible for these devices 
not being able to achieve Fano- limited spectral resolution 
in the optical and soft X-ray band 2-4 . 

Traditional X-ray microcalorimeter geometries of- 
ten employ thin-film membranes 5 or narrow silicon legs 6 
to isolate the temperature sensor from the substrate. 
In these geometries athermal phonon energy loss is of- 
ten substantially reduced or almost eliminated. With 
a membrane geometry, athermal phonons are preferen- 
tially thermalized or reabsorbed rather than lost from 
the detector due to the physical geometry. An energy 
resolving power of over 3000 has been achieved while de- 
tecting 6 keV photons 7 , and over 4000 while detecting 
100 keV photons 8 . While this very high spectroscopic 
performance is achievable with this type of device geom- 
etry, the low thermal conductance of the membranes can 
lead to detectors that return to the temperature of the 
thermal bath with a relatively slow time constant, lead- 
ing to lower count-rate capability than desired. There 
are many applications where high count-rate capability 
and high spectral resolution is desired 19 . Another draw- 
back to using membrane supports is that large arrays of 
fine-pitch and close-packed pixels become more difficult 
to design and fabricate. 

To produce a fast microcalorimeter that does not 
have degraded energy resolution a geometry was fabri- 
cated with a TES sensor on a solid substrate which is 
attached to an X-ray absorber supported by attachments 
points (“stems”), that supports it above the sensor and 
substrate. By making these contact points small, the ef- 
fects of athermal energy loss were minimized 12 . It was 
found that for this geometry, the energy resolution de- 
graded as the area of contact of the “stem-contact” in- 
creased. The detector geometry with the smallest contact 
area between absorber and solid substrate was not only 
fast, but also demonstrated the best energy resolution 
seen by any low-temperature detector at 6 keV of 1.6 eV 
[FWHM] 12 . In similar devices where the absorber was 
twice as thick, the dependence of energy resolution on 
contact area was not apparent. Thus there was strong 
evidence for the energy resolution of this type of detec- 
tor being highly dependent on the amount of athermal 
phonon loss, which is very geometry-dependent. It is 
the understanding of the important physical processes 
that takes place when athermal phonons are generated 


in metal films that motivates this work. With this un- 
derstanding, it will be possible to predict the effects on 
performance for a wide variety of absorber and stem di- 
mensions. It will also be possible to optimize the geome- 
try and material properties according to the requirements 
of the different detector types, taking into account the 
need to simplify fabrication processes as much possible. 

In Section II we describe two experimental results 
that are affected by athermal phonon loss. In the first, we 
consider the performance for when X-ray are stopped in a 
microcalorimeter sensor film on a solid substrate without 
an absorber. In this type of geometry, the performance 
is known to be poor for decades due to athermal phonon 
loss. However we use the width and shape of the spectral 
response from the almost mono-energetic source as an 
important measurement for comparison with our model. 
In the second we look at signals from a microcalorime- 
ter with a similar film on a solid substrate with an ab- 
sorber attached. In this latter case, by comparing the 
signal profiles when X-rays are stopped in the absorber 
to those that are stopped in the sensor, we are able to 
measure important parameters in the energy transport 
process. For both of these experiments the sensor is 
a magnetic penetration-depth thermometer (MPT) 11 , a 
relatively new type of thin-filni thermometer. It is partic- 
ularly useful for these studies because of important sen- 
sor properties that we describe in Section II that ensure 
the energy loss observed is dominated by the processes 
modeled. In Sections III-V we describe a framework for 
studying how athermal phonon loss can affect the perfor- 
mance and characterization of different detector designs 
on solid substrates. 

Prior to this paper there have been no direct mea- 
surements or calculations of energy loss and fluctuations 
in this loss that will influence measured spectral line- 
shapes. For the high resolving power of TES sensors fab- 
ricated on suspended membranes this was not important. 
For STJs the energy loss can only be evaluated indirectly 
through line shifts which must be correlated with certain 
sharp features in the dependence of photon absorption 
coefficient on energy either due to peaks associated with 
onset of inter-band transitions for optical photon range 
or in X-ray absorption edge spectroscopy for soft X-rays. 
In both cases the abrupt change in photon absorption 
coefficient and distribution of absorption sites result in 
strong variations in energy loss and corresponding line- 
shifts. However, indirect measurements cannot yield suf- 
ficiently accurate results because of practical experimen- 
tal details such as the poor quantum efficiency of thin 
metal films, and often complicated fine X-ray absorption 
spectra. The direct measurements of energy loss from 
thin layer systems through athermal phonons therefore 
is of great importance. The results are relevant to a wide 
variety of thin-film metallic sensor geometries, where the 
sensors can be transition-edge sensors (TES) 9 , metallic 
magnetic calorimeters (MMC) 10 , magnetic penetration 
thermometers 11 , MKIDs, STJs, and nano-bolometers. 
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FIG. 1: (Color online) Schematic depicting the components 
of an MPT microcalorimeter (left) and an assembled pixel 
(right). 

II. EXPERIMENTAL DETAILS 

In this section we describe two experimental results 
that provide information about the energy loss that oc- 
curs when an X-ray is stopped in a thin metallic film 
on a solid dielectric substrate, prior to the film reach- 
ing quasi-thermal equilibrium. For both of these ex- 
periments, the quasi-thermal signal is measured using a 
magnetic penetration thermometer (MPT) 11,13,14 . This 
type of thermometer utilizes the temperature-dependent 
diamagnetic susceptibility of a superconducting film. A 
schematic diagram of the two experiments is depicted 
in Fig. 1, which shows the layout of an MPT-based X- 
ray microcalorimeter. When an X-ray is stopped in the 
gold absorber the entire absorber and the MPT super- 
conducting film quickly come into thermal equilibrium. 
The MPT-film is a bilayer consisting of molybdenum and 
gold. The temperature of the microcalorimeter ideally 
rises according to 

AT = AE/ Ctotai (1) 

where A E is the energy of the X-ray, and C to tai is the 
heat capacity of the absorber plus the MPT bilayer. The 
temporal profile of a typical event where an X-ray is 
stopped in the absorber is shown in the inset of Fig. 2. 
After the first ~ 1 /is the electron and phonon systems 
are in equilibrium, but the temperature of this system is 
changing. Thus we consider the device to be in quasi- 
thermal equilibrium. After the temperature rises the 
microcalorimeter decays back to the temperature of the 
heat bath with an exponential decay time constant given 
by Ctotai /G, where G is the thermal conductance of the 
MPT film to the substrate. G is determined by the 
Kapitza conductance of the MPT film to the dielectric 
substrate below. The signal is read out with a pick-up 
coil that is beneath the MPT superconducting film. This 
pick-up is connected to the input coil of a SQUID forming 
a purely superconducting loop. As the temperature rises, 
the MPT diamagnetic susceptibility is reduced, and the 
inductance of the pick-up coil increases. The inductance 
( L p (T )) is sensed from current seen through the SQUID 
input (A) according to 

j _ -h [L P { oo) - L p (T)} 

1 L P (T ) + Li 



Energy [eV] 

FIG. 2: (Color online) Spectrum from MnK a X-rays de- 

posited directly into MPT bilayer. The natural line shape 
of the MnK X-rays is shown as the narrow solid lines un- 
derneath the data points. Inset is the temporal profile of a 
typical MnK a pulse. 



FIG. 3: (Color online) Measured SQUID signal as a function 
of temperature for MPT without absorber attached. Ib=l-0 
mA, producing an anisotropic magnetic field in the region of 
the MPT bilayer. 

where Li is the sum of the SQUID input inductance and 
stray inductance, and is the near constant bias cur- 
rent supplied to the parallel combination of the SQUID 
input coil and the MPT pick-up coil, provided by a large 
ballast inductor in the bias circuit. Ib is actually a per- 
sistent dc current that is applied with a procedure that 
ensures that almost all of the current flows through the 
pick-up coil 15 . The magnetic field produced by the pick- 
up is non-isotropic and has been investigated in detail 
by Fleisclimann et al. 10 . In Fig. 3 we see an example 
of a SQUID signal as a function of temperature with an 
anisotropic field, which is approximately the temperature 
dependence of the MPT diamagnetic susceptibility. Un- 
like many other microcalorimeters, the MPT film has no 
electrical connections, and so the possibility of electric 
connections that could allow energy loss through these 
contacts is eliminated 16 . 
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These MPT microcalorimeters are fabricated using 
standard microfabrication techniques 17 . There is a 
250 nm thick AI2O3 layer that provides isolation between 
the Si substrate and the first Nb layer, which forms the 
pickup coil. To isolate the MoAu sensor from the Nb 
pick-up coil, we anodize the Nb forming a 50 nm thick 
Nb2C>5 layer and then deposit a 250 nm thick layer of 
AI2O3. A second Nb layer is then deposited to form low- 
inductance microstrip leads in all regions of the circuit 
except the pick-up coil. Deposited over the area of the in- 
sulated pick-up coil is the MoAu superconducting bilayer. 
The lower layer of the bilayer is 50 nm thick Mo, and the 
upper layer is 330 nm thick Au. In the first and second 
experiments the pick-up coil had different geometries and 
calculated inductances of 0.76 nH and 0.53 nH respec- 
tively. The MPT chip is connected to a separate SQUID 
chip with aluminum wirebonds, which contribute a stray 
inductance, typically measured to be «1.6 nH. Low tem- 
perature characterization of our devices is performed in 
an ADR capable of reaching a minimum temperature 
of 30 rnK. Our readout is accomplished using a two- 
stage SQUID provided by PTB-Berlin 18 . The first stage 
SQUID has an input inductance of 1.8 nH, an input cou- 
pling of 5.4 /tA/$ 0 , and the overall readout has a broad- 
band noise level that is less than 0.4 /t^o/a/Hz at our 
operating temperature referred to the first stage SQUID. 
The SQUIDs are read out using Magnicon SQUID elec- 
tronics. Further details of our experimental procedure 
can be found in Rotzinger et al 15 . 

In the first experiment, where there was no X-ray ab- 
sorber, the X-rays were stopped directly in the MPT film. 
The main part of the X-ray spectrum of pulse sizes is 
shown in Fig. 2. These signals are from 5.9 keV Mn K a 
X-rays that are produced when a radioactive 55 Fe source 
undergoes electron capture; signals from Mn were also 
observable at a lower count rate. The ratio of Mn K a and 
K /3 signal sizes, 500 eV apart, was used to calibrate the 
responsivity. The Mn K Q X-rays are practically mono- 
energetic on the scale of the signal-size variation seen 
in this measurement. Based upon the signal and noise, 
the detector resolution for this measurement was 4.1 eV 
[FWHM] and the Mn K a complex extends only over a 
range of 10 eV, as seen in Fig. 2. We are able to compare 
the distribution of pulse sizes with those we would ex- 
pect from mechanisms that we can model to learn about 
energy loss processes. Fig. 3 shows the SQUID output as 
a function of temperature for this particular device. For 
this measurement with /{,=1.0 mA, it was necessary to 
operate at a temperature of 115 rnK. This ensured a re- 
sponsivity that was only slightly non-linear with a signal 
gain that could be accurately calibrated. It also ensured 
that the response from the X-ray signal was not too large 
to maintain the flux-locked loop of the SQUID read-out 
from the ~ 12 mK temperature rise that is seen on the 
time-scale of the signal rise-time which is ~ 0.2 [is. 

In the second experiment, with a different geometry 
MPT sensor coil, the X-ray absorber is attached and the 
majority of X-rays are stopped in this absorber. The 


absorber is a thin gold foil that is suspended above 
and coupled to the MoAu film; the absorber dimen- 
sions are 250x250x2.8 /1m 3 . In Fig. 1 the connec- 
tion points are shown schematically as solid “steins” ; 
for these devices the contact “steins” were small regions 
where the absorber material is lowered to contact the 
sensor and substrate with small contact points that are 
3.5 /im x 3.5 /im. There was a single contact point onto 
MoAu, and four contact points onto just the dielectric 
substrate that were introduced to add mechanical sta- 
bility. The total area of contact of these stems is just 
0.1 % of the area of the MoAu. The superconducting 
transition temperature was lower for this experiment and 
the temperature of operation of 35 mK is at the steep- 
est part of the MPT transition curve for 1^=0. 38 mA 13 . 
When X-rays are stopped in the absorber, i.e. when the 
whole energy of the X-ray was retained in the spectrom- 
eter following its absorption, the measured spectral res- 
olution of MnK Q complex was 4.3 eV [FWHM] 13 . This 
spectral resolution level is consistent with the resolution 
expected based upon the signal and noise. It is only 
limited by the microcalorimeter thermodynamic fluctu- 
ations and the SQUID amplifier noise level, given the 
temperature responsivity of this device under the oper- 
ating conditions. Thus for this design, with a relatively 
small contact area between the absorber and solid sub- 
strate, the energy resolution is apparently not affected by 
the de-grading effects of athermal phonon loss. In this 
respect it is close to “ideal” . 

A small fraction of the X-rays passed through the ab- 
sorber and were stopped in the MPT Him. Based upon 
the geometry of the device and the known attenuation 
lengths of the different materials we estimate that 1.55% 
of the 5.9 keV X-rays incident on the absorber actually 
pass through it and are stopped in the MoAu sensor. 
The measured fraction of 1.85% ± 0.23% is within the 
range of what we would expect given the known accu- 
racy of absorber thickness (~10%). The events where 
the X-rays are stopped in the bilayer have a very differ- 
ent temporal profile, as shown in Fig. 4. We can see a 
very fast signal rise (< 1 fi s) followed by an almost as fast 
initial decay (~ 2.7 [i s) as the MPT film comes into quasi- 
thermal equilibrium with the absorber, and then the sig- 
nal decays with the slower time constant associated with 
the microcalorimeter sensor plus absorber heat capacity. 
What is interesting is that although the slow time decay 
constant is identical for X-ray events into the absorber 
and the MPT film, the signal is measurably smaller when 
absorbed in the MPT film. By fitting to the regions of 
the temperature decays where there is just a single expo- 
nential decay we determined that the MPT absorption 
events are 3.8 ± 0.1% smaller. It is this difference in 
signal size that we interpret as mostly being the energy 
that is lost to the substrate through athermal phonons. 
One alternative mechanism for losing heat is due to the 
additional thermal conductance of heat out of the mi- 
crocalorimeter due to the MPT film being at a higher 
temperature for a short duration. Knowing the heat ca- 
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FIG. 4: (Color online) Temporal profiles of the MnK a X-ray 
signals stopped in the absorber and in the bilayer. Inset is 
a plot of the same signals on an expanded time-scale, with a 
logarithmic scale for the signals. 

pacifies of the MPT film and the absorber, and knowing 
the thermal conductances, we have modeled this heat 
loss path and determined that the additional amount of 
heat that is conducted away thermally is approximately 
0.4±0.1%. Thus we attribute the remaining 3.4% of the 
difference in signal heights as being due to energy loss via 
athermal phonons. In addition to the obvious difference 
in pulse shape over the first 2 /zs there is a small depar- 
ture of the two curves from simple exponential decays 
for approximately the first 100 /zs after the absorber and 
bilayer have come into thermal equilibrium (represented 
by the dashed lines in Fig. 4). 

III. THEORETICAL SIMULATION OF ENERGY 
LOSS AND SHAPE OF THE SPECTRAL LINES 


The use of the term “athermal” phonons refers to a 
situation where phonons as well as all other excitations 
are far away from thermal equilibrium. This term how- 
ever is too broad to properly describe the important fea- 
tures of a non-equilibrium phonon distribution. Nonethe- 
less we adopt this terminology while discussing general 
concepts. This is not misleading because in all situa- 
tions we characterize interacting phonons and electrons 
participating in the energy down-conversion cascade in a 
complete way using time evolving solutions of the appro- 
priate kinetic equations. 

The energy down conversion following the absorp- 
tion of an energetic photon in a metal can be viewed 
as evolving through three distinct stages 2 * 4,20 . Stage 
I starts from the ejection of a photoelectron and pro- 
ceeds through electron-electron (e-e) interactions, finish- 
ing with the formation of a highly non-equilibrium dis- 
tribution of electrons and holes in the conduction band 


occupying a tiny volume in the vicinity of an absorption 
site. During the stage I the whole of photon energy is 
transferred from photoelectron to the ensemble of sec- 
ondary, ternary, etc. electronic excitations with rapidly 
diminishing mean energy. 

Stage II takes over when the mean energy of 
electronic excitations becomes sufficiently low for the 
electron-phonon (e-ph) interaction to determine further 
momentum and energy relaxation processes. This hap- 
pens below a threshold E jf, i.e. e < El where e is the 
electron energy. For typical metals El far exceeds the 
Debye energy, El tip. During this stage electrons and 
holes mostly relax via emission of phonons. As a result 
a bigger fraction of photon energy is transformed into 
the energy of the evolving phonon distribution. Stage 
II is completed once the mean excitation energy reaches 
the low energy threshold S2i which is determined from 
the relation T e - p h{Q i) = T p h- e (^i)- Below this thresh- 
old (stage III) a rapid energy conversion into predomi- 
nantly electronic excitations takes place. The scenario 
described is correct provided the electron-electron inter- 
action in the whole region e > Dj remains weaker than 
the electron-phonon interaction with spontaneous emis- 
sion of a phonon, r ee (e) > r e _ p h(D i). Whether this 
is true depends upon material parameters and on the 
amount of disorder, which has a strong effect on the 
electron-electron interaction. We will consider both cases 
keeping the notation fli for the end point of the stage II 
accounting possible role of electron-electron interactions. 

During stages I and III practically all the initial pho- 
ton energy is stored within the electronic system and 
cannot be lost from the detector. In contrast with elec- 
tronic excitations, phonons can easily escape through 
metal/dielectric interfaces. Therefore the only stage dur- 
ing which a significant energy loss may occur is stage II 
where athermal phonons carry most of the photon energy. 

The simplest planar geometry that we will consider 
is shown in Fig. 5. The only possible phonon escape 
interface is at z = —d/2. The him is supposed to be 
grown on a solid substrate occupying the half space z < 
—d/2. Since z = +d/2 is a vacuum interface this is not an 
escape interface. There are several phenomena associated 
with escape of athermal phonons 3 . 

First, the energy loss, Ei oss (z), is a function of the 
coordinate of an absorption site relative to the escape 
interface. We only show the z-coordinate dependence, 
where the z-axis in Fig. 5 is in the direction normal to 
the him. Lateral coordinates are not important as long as 
photon absorption does not occur very close to the edges 
of the him. However, the z-dependence is important de- 
termining the distance to the escape interface relative to 
mean free paths of athermal phonons. 

Second, only phonons propagating inside the critical 
cone can pass through the interface, otherwise they are 
internally rehected. The existence of this restriction on 
transmission introduces statistical fluctuations to the loss 
of energy due to fluctuations in number of phonons emit- 
ted into the critical cone. The statistical fluctuation in 
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FIG. 5: Sensor film geometry 


the number of escaping phonons contributes to a Gaus- 
sian fluctuation in the energy loss, and can be viewed 
as an specific “athermal phonon noise” , and thus can be 
added in quadrature with other statistically independent 
noise sources. These fluctuations are important for opti- 
cal and near infrared photons, but their effect becomes 
progressively small for X-rays 20 . 

Finally, since the energy loss depends on the coordi- 
nate 2 the distribution P(z) of the X-ray absorption sites 
is important factor determining the spectral line shape. 
The line shape depends on the X-ray energy and the an- 
gle of incidence. 


A. Energy loss due to phonons of the first 
generation 

The phonon distribution of the first generation is 
established during the very fast electron down-conversion 
Ei — > flp and was found in explicit form in our earlier 
work. Energy loss due to the escape of phonons of the 
first generation created at the point of absorption at Zq 
has the form 20 


E, „»»,,(*,) = 4E g; 1/2 + *,/<*) x 


m = 0 


1 H" Sq 


m, 0 


de 


y„ Tr D {n~ D 

[1 - exp (irrnr - d/l ph - e (e)£,)] 
l + m 2 TT 2 l 2 ph _ e (e)^/d 2 


lph—e (c) 


( 3 ) 


where we denote k(x) = exp(— x) sinh(a;)/x, and define 
£ 2 = 7r 2 Dtdc/2d 2 , D is the diffusion coefficient of down- 
converting electrons (e > f Id), tdc is the duration of 
E* — i Qd stage and l p h- e {t) is the phonon mean free 
path with respect to electronic re- absorption. The trans- 
mission coefficient r](£) is a function of incidence angle, 
and £ c = cosd c , where i9 c is the critical angle. 

We simplify the general expression (3) further. In 
what follows we use the linear dispersion relations for 
phonons taking wq = cq, where c is the velocity of sound 
and q is the phonon quasi-momentum. In many low tem- 
perature detectors such as TESs, MPTs and MMCs the 
sensor film is either a normal metal or a superconduc- 
tor operated close to the critical temperature T c . Even 
in most of other detectors, such as STJs and MKIDs, 
which are operated at T < T c the lower energy thresh- 
old of the phonon stage of down-conversion cascade sat- 
isfies the condition Q-| A. Under this condition both 


electronic excitations and the matrix elements for the 
electron-phonon interaction are essentially the same as in 
a normal metal. Therefore we will consider the asymp- 
totical limit A — > 0. In this limit the expression for the 
phonon lifetime due to re-absorption by electron system 
in metal is 


_i 47r./V(0)a 2 (e)e 

W e) = HN 


( 4 ) 


where a 2 (e) is a coupling strength for electron-phonon 
interaction entering the Eliasliberg function a 2 (e)F(e), 
F(e) is the phonon density of states, iV(0) is the density 
of states per spin at the Fermi level and N is the den- 
sity of atoms. Although a 2 (e) is a function of phonon 
energy and may exhibit variation over the phonon spec- 
trum, the model becomes especially simple if we can use 
for it some averaged value. We thus introduce another 
model assumption taking a 2 (e) as a constant, so that 
T ph- e (t ) ex e. In this case the last integral in the expres- 
sion (3) has the form 


Y{m,P,£) 


p 4 e 



dxx^ 


1 — exp(i?n7r — x) 
X 2 + TO 2 7T 2 
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where /3 = l p h- e {^D)/d. This integral can be explicitly 
calculated. It has the simplest form in the asymptotic 
limit /? — ► 0. Neglecting small terms of the order of 
exp(— 1//3) <lwe obtain 


Y(m,0,£) = t {l — 3?n 2 7r 2 /3 2 £ 2 — 3/3 3 £ 3 (— l) m (2 — m 2 n 2 ) 

O 


+3m 3 7r 3 /3 3 £ 3 


arctan 



+ Si(m7r) 
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where Si(;r) is the integral sine function. Integration over 
£ in the expression (3) is within the limits between £ c 
and 1. For acoustically soft metals on rigid substrates 
£ c is close to unity and because the function F(?n,/3,£) 
depends smoothly on £ it can be evaluated at £ = 1. 
Introducing fj = f* d££r/(£)/ J ^ d££ and p = sin 2 (0 C / 2) 
we finally obtain the energy loss due to escaping phonons 
of the first generation in the form 


Ei oss,l (z 0 ) = dEpp(l +£ c ) x 
^k(to 2 C 2 )costott 0 + ^ Y(m,P, 1) (7) 

In expressions (3) and (7) there are two separate 
factors that determine the magnitude of loss: the proba- 
bility of phonons to reach the escape interface (depending 
on intrinsic material properties); and transmission char- 
acteristics of the interface. 


B. Energy loss during Qd — ¥ fii sub-stage 

The energy loss from the first generation of phonons 
may not be the dominant mechanism for energy loss be- 
cause their mean energy is very high and so their mean 
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free paths are relatively short. For the same reason the 
first generation of phonons rapidly decays and after one 
- two generations asymptotically converges to the auto- 
model solution of the kinetic equations describing the 
interaction between electrons and phonons. The details 
of this solution were discussed in our earlier work 21 . 

Each subsequent phonon generation is characterized 
by a lower mean energy. In the first phonon generation 
the majority of phonons occupy states in the region of 
the highest density of states, i.e. close to the Debye en- 
ergy. Thus the lifetime of first generation of phonons 
is roughly T p h- e (^D)- The absorption of one of these 
phonons creates two electronic excitations with mean en- 
ergy of ^ljj/2. The lifetime of electrons in metals with 
respect to the spontaneous emission of a phonon is 


strong inequality fli is not fulfilled that contribu- 

tions from second phonon generations may be needed. 
An expression for Ei oss q d ^ > .q 1 can be written as 





12Z(0) 
11Z(0) + 3 


Pa 



(9) 


where 

r 1 dz 

g(x) =X —f(z) {Ei(l , z{x - 1)) - Ei(l, 1 - z) 
Jl/x z 

+ exp(z) [Ei(l, 1) - Ei(l, xz)] + In ^ j , 

Ei (a, z) is the exponential integral and 


-i 


' e-ph ( e ) — 


27 r 

nz{ o) 


dfla 2 (fl)F(fl) (8) f(z) = 1-z - — jz 2 cos \[2 In z — 7\Pl sin V2 In 2 — l| 


where Z{ 0) is the renormalization factor of a metal. For 
e > Qd we obtain T~ x ph (e) = T~} ph (Vl D ) = t s and this 
time does not depend on electron energy. For e < Hd and 
assuming a 2 is a constant value, r~ - P h ( e ) « e3 for the 
Debye model of phonon spectrum. From this estimate it 
follows that the intermediate electrons (holes) will diffuse 
during their average survival time, i.e. T e - p h{^D/ 2 ) = 
2 3 t s . 

Another important factor is the evolution of the elec- 
tron diffusivity. During the E* — v Qd stage the electron 
momentum relaxation is determined by the combined ac- 
tion of elastic scattering and electron-phonon scattering 
with the rate r _1 = t~ x + rj 1 . The intermediate elec- 
trons have a typical energy of they emit lower 

frequency phonons of the second generation and there- 
fore are scattered less efficiently. For these electrons 
T _ 1 ( fl £)/ 2 ) = t~ x + t~ 1 / 2 3 so that D' > D, where D' 
is their diffusion coefficient. This diffusion coefficient is 
faster than for electrons with e > Qd resulting in a much 
greater spread of phonons of the second generation. The 
extent of spatial spread can be evaluated by comparing 
the factor Dtd c for the first generation of phonons with 
8 D't s for the second. 

The accurate evaluation of the contribution of the 
second generation phonons to energy loss, fluctuations 
and spectral line shape is difficult. As we will see below, 
most of the energy loss occurs at an advanced stage of the 
fin — > cascade. Starting with the decay of phonons of 

the first generation, the phonon distribution rapidly con- 
verges to an asymptotic automodel solution 2 . The auto- 
model solution at least partially takes into account con- 
tributions from the second phonon generation. The dom- 
inant role of fin — ► Hi stage in energy loss arises because 
as the mean phonon energy decreases, more phonons are 
likely to reach the escape interface. In contrast, phonons 
of the first generation make dominant contribution to 
fluctuations of the energy loss and asymmetry of the spec- 
tral line shape. In metals with a long fljy — >• fli stage, 
i.e. flu A> lb , we expect the automodel solution to be a 
good approximation. It is only in the situation when the 


In equation (9) we kept only spatially homogeneous con- 
tributions and neglected all spatially varying terms in the 
general solution 21 , which are not important. The order of 
magnitude of the neglected terms is d 2 /Tr 2 DT p h- e (^D), 
which is much less than unity for all realistic situations. 

Combining all results we finally arrive in the total 
energy loss from our sensor him 

Ei OSS (z 0 ) = Ei oss,l (z 0 ) + Ei 

OSS,Q.D—*Ql ( 10 ) 

where the first contribution is strongly dependent upon 
the z-coordinate of absorption site, while the second con- 
tribution is highly insensitive to such a dependence. The 
magnitudes of both contributions are determined by the 
effective transmission coefficients for the first and sub- 
sequent generations of phonons, and the extent of the 
Qd — > fli sub-stage of the stage II down-conversion cas- 
cade. In an ideal him the second contribution can be 
the largest and is determined by the intrinsic proper- 
ties of the him material. In real films the duration of 
fin — > Di sub-stage may become signihcantly shortened 
by a strong elastic scattering enhancing electron-electron 
interaction. The dependence of £) OSS) n D _>.Qi on the pa- 
rameter tti/Qc, is illustrated in Fig. 6. This dependence 
is strong, so that the relative contributions of the first 
and subsequent phonon generations to the total energy 
lost is sensitive to the strength of the electron-electron 
interaction. 

Finally the integrands in expressions (3) for 
Eioss, 1 ( 20 ) and the function g(x) for Ei oss ^q d ^q 1 in (9) 
give the spectral distribution of athermal phonons escap- 
ing from a metal him into a substrate. There are many 
physical situations where this information may be im- 
portant. Here we mention just two examples. The first 
is in calculating the amount of energy back-how into a 
thermalizing metal him. Due to the strong elastic scat- 
tering in the substrate (either isotope-related or due to 
disorder), the back-scattered phonon hux depends upon 
the energy of athermal phonons. It can be evaluated pro- 
vided that the spectral distribution of escaping phonons 






FIG. 6: (Color online) The dependence of Ei OSSi sn D - > n 1 on 
low threshold fii /f2o 


is known. In special cases this energy back-flow may ul- 
timately affect the energy loss from the film. A second 
example is in quantifying the thermalization or energy 
down-converting properties of thin metal films grown on 
the surfaces of semiconductors or dielectrics where there 
is an incident flux of high energy phonons. In this situa- 
tion the excitations that are generated within the metal 
film are induced by the incident athermal phonon flux. 
With a slight modification to the initial phonon distri- 
bution in the film the quantitative relations for spectral 
transformations of the reflected flux relative to the inci- 
dent flux can be derived. 


C. Spectral line-shape 

In the previous subsections we presented several im- 
portant asymptotical results derived from the general 
expressions as preliminary material for the original de- 
velopment below. In this subsection we start discussing 
modeling of the spectral line-shape. 

This topic in the area of low temperature detectors 
so far has been restricted to very limited number of cases, 
namely, detectors with a spatially-homogeneous response 
where the spectral line-shape is Gaussian, and spatially 
inhomogeneous STJs where the standard approach was 
fitting the experimental data with a skewed Gaussian al- 
lowing the variance to be dependent on energy. System- 
atic attempts to study the role of spatial inhomogeneities 
in lateral directions and simulate the effects of quasiparti- 
cle losses at the edges or in the leads or local defects have 
been reported in several works both experimental 22-24 
and theoretical 25-28 . The only theoretical simulation of 
the spectral line-shape 28 in laterally inhomogeneous de- 
tectors has revealed that spectral line-shape is sensitive 
to the details of energy loss. 

In the nearly ideal detectors which we discuss in 
this paper we focus our attention on the fundamental 


^-dependence of Ei oss and ignore the less fundamental 
potential effects of lateral inhomogeneity. The only dis- 
cussion of the role of “vertical” inhomogeneities on wors- 
ening of the FWHM, when fitting the experimental lines 
with a skewed Gaussian, is available in Kozorezov et al., 
20 07 20 . We will demonstrate below in Sections IV and 
V that further development of the model aimed at sim- 
ulation of spectral line-shape of nearly ideal detector is 
possible and desirable. It becomes very informative be- 
cause of the sensitivity of detector response to details of 
energy deposition and loss. 

In the ideal situation when there is no energy loss, 
Eioss = 0, the normalised spectral lineshape S(E) is 
Gaussian, J^° oo dES(E) = 1, 


S(E) 


, — exp 
\Z2na 


(E — £’o) 2 l 
2<r 2 


( 11 ) 


The line is centered around photon energy Eq, and it’s 
full width at half maximum is 2y2 In 2er. If Ei oss {zq) ^ 0, 
then the lineshape is 


S(E,z 0 ) 


, — exp 
\Z2na 


(E — Eq + Eioss(zo )) 2 
2cr 2 


( 12 ) 


provided that an absorption takes place at a site with 
fixed coordinate zq. If P{zq,Eq) is the normalised 
(f dzoP(zo, Eq) = 1 ) distribution function of absorp- 
tion sites at a depth zq f° r a photon energy Eq then the 
spectral lineshape is given by a convolution 


S(E) = f dz 0 P(z 0 ,E 0 )S(E,z 0 ) 

J detector 

The lineshape due to absorption events in the sensor film 
is 


S(E) = 


fd/2 


/ -d/2 


dz 0 P s (z 0 ,E 0 )S(E,z 0 ) 


(13) 


where the distribution of absorption sites in the sensor 
for X-ray incident onto z = +d / 2 interface is 


Ps(z 0 ,E 0 ) 


exp[(z 0 - d/2)/L(E 0 )\ 
L(E 0 )[1 - exp(— d/ L{Eq))] 


(14) 


where L(E 0 ) is the 1/e photon absorption depth. Finally, 
if the source is not monochromatic an extra straightfor- 
ward convolution with its spectral density is needed. 


IV. MODELING SPECTRAL LINE-SHAPE AND 
ENERGY LOSS FROM X-RAY ABSORPTION 
EVENTS IN MO/AU FILM 


In this section we concentrate on modeling the spec- 
tral line-shape and energy loss from X-ray absorption 
events in a Mo/Au film. In soft X-ray range the role 
of fluctuations on spectral line-shape was shown to be 
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generally small 20 . For soft X-rays at 6keV using the ex- 
pression for fluctuation factor 20 , J, we estimate the vari- 
ance due to fluctuations in number of escaping phonons 
of the first generation for Mo/Au sensor to be a f <1.0eV. 
For all purposes of this work this is small and can be ne- 
glected. 

The two different experiments described in Section 
II allow us to probe different physical parameters of the 
film. In particular, energy loss is mostly sensitive to the 
duration of stage II and transmissivity of the escape in- 
terface, but not to the size and position of the initial 
phonon bubble (first generation of phonons) relative to 
the escape interface. In contrast, the spectral line-shape 
is sensitive to both the size and the position of a phonon 
bubble. Another aspect to consider is that phonon trans- 
mission may be strongly assisted by scattering if there is 
significant roughness of the interface. Fitting the model 
predictions to both sets of data potentially reveals details 
of the energy that is transferred across the escape inter- 
face by athermal phonons over a broad frequency range. 
It also allows us to evaluate the strength of the electron- 
electron interaction relative to the electron-phonon inter- 
action in metals, both in high excitation energy range, 
from few eV down to the Debye energy above the Fermi 
level and at much lower energies close to equilibration. 
These insights, apart from their obvious practical impli- 
cations for single photon low temperature detectors, also 
have more fundamental importance for characterizing the 
interactions in highly excited metals and the equilibra- 
tion process over a much wider spectral range than those 
previously accessible with other experiments. 


A. Definitions of down-conversion stages I and II 
in Au and Mo 

In our modeling, we start from the time during stage 
I when the e-ph interaction becomes stronger than the e-e 
scattering. This happens at the threshold Ei determined 
from the condition r~_ e (Ei) = r” 1 . In our earlier work 2 
we considered the metal to be sufficiently pure that for E\ 
in the eV range above the Fermi energy the e-e scattering 
could be considered as in ideal metal and not affected by 
disorder. However, in thin films with residual resistivity 
ratios (RRR) less than 10 it is often the case that elastic 
scattering times for electrons, Ti mp , either due to bulk 
scattering or diffuse scattering at surfaces or interfaces, 
are of the same order of magnitude as r s . Correspond- 
ingly, at the energy Ei , the electron-electron scattering is 
significantly affected by elastic scattering. In this energy 
range single electron-phonon scattering events are quasi- 
elastic and also contribute to elastic scattering. The ef- 
fect of elastic scattering on electron-electron interactions 
depends on the system dimensionality 29 . The criterion 
for the normal direction of the film of thickness d to not 
contribute to the dimensionality 30 is d < \/TiD/e. For 
the lowest energies of interest, e > Qi , the system be- 
comes 2-dimensional in the context of e-e scattering only 


if d < y/fiD/ fii. Thus we expect that even for small 
thresholds <C II d unless the film is very thin it is 
still 3-dimensional for most of the f Ip -* fli interval 
so that T“ 1 (e) oc e 3 / 2 . At the same time both impu- 
rity and electron-phonon scattering times in this energy 
range may be taken to be independent of energy. As a 
rough estimate we scale r (e) up in energy according to 
the e 3 / 2 law until it becomes as big as t// + r” 1 . We 

thus re-define E\ so that in less pure metals it sets the 
threshold below which electron momentum scattering is 
controlled by impurities and phonons 


Ts "b Emp 1~e— e(El) e, ideal {El ) 


= a 


El 

TiEp 



15) 


where Ep is the Fermi energy, and r s is the radius of a 
sphere in atomic units which encloses one electron charge, 
and a > 1 is a numerical factor that takes into account 
that at Ei the disorder scattering slightly enhances the 
e-e scattering relative to the ideal situation. The last 
term in (15) is the expression for the e-e scattering rate 
in an ideal metal 31 . 

At the end of stage I at an energy E J the electron- 
phonon scattering also becomes a dominant source of 
energy relaxation. The latter can be found from the 
equation 32 


Tee(El) = ||t s (16) 

The meaning of this balance is that over the duration of 
the time equal to T ee (El), the probability of having the 
electron energy relaxed by a factor 3 in the e-e scatter- 
ing process (splitting the initial energy El equally be- 
tween the three electronic excitations after e-e intera- 
tion) is the same as for the process of sequential emis- 
sion of the 2El/3flp Debye energy phonons. Taking 
r” 1 = t~J-{Ei) (El/Ei) 3 ^ 2 and using (15)-(16) we ob- 
tain 


e* = n D 


13.67 TiEp \ 3/1 ° 

K ) 


rV 2 T 


1 + t s /t i 


imp 


1/10 

(17) 


To estimate the duration of the down-conversion 
stage from E\ down to flp we first evaluate the du- 
ration of Ei El, i.e. tdc{Ei El). Assuming 
that with each e-e scattering event the energy of elec- 
trons decreases by a factor of 3, we must sum the ge- 
ometric progression with the scale factor 3 3 / 2 . Thus 
t-dc.E^Ei < T ee {El)/{ 3 3/2 - 1). This is an upper limit 
to the estimate because the geometric progression must 
be truncated due to finite ratio Ei/El and weaker par- 
allel phonon emission also takes place. Below El the 
relaxation down to a level El / 3 occurs over the duration 
T ee (El)/2 with the e-e and e-ph scattering processes at 
El contributing equally. Further relaxation from El /3 to 
fin lasts for another T ee (El)/ 2 only taking into account 
phonon scattering. Combining all of these estimates we 
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have t dc < (l/(3\/3— 1) + l) r ee (El) » 1.24 T ee {E{) » 

(E$/n D )r s . 

More detailed studies of the down-conversion stages, 
I and II, require estimates of the electron-phonon cou- 
pling strength in the bilayer. In practice this is very dif- 
ficult to analyze. Indeed, gold belongs to a class of weakly 
coupling metals, while molybdenum - to the intermedi- 
ate coupling metals. We will make estimates for Au and 
Mo separately before hypothesizing about the appropri- 
ate parameter values for the Au(330nm)/Mo(50)nm bi- 
layer. 

We first estimate t s . Taking t s = T e - p h(^D) from 
(8) and using the experimentally measured Eliasherg 
function a 2 (H)E(H) for gold 33 we obtain r s « 0.13ps. A 
similar estimate for molybdenum based on the calculated 
Eliashberg function 34 yields t s s=s 0.030ps. Thus, td c in 
Mo is likely to be a factor of 4 faster than in Au assum- 
ing roughly the same ratio E^/tlp), (in Mo f Ip = 33meV, 
and combination t s £1 2 d in (17) is nearly the same for both 
materials). 

In order to evaluate T p h- e ( II) we use expression (4). 
Assuming a 2 (H) to be a constant value and expressing 
it in terms of the first reciprocal moment of Eliashberg 
function A (mass enhancement factor), in the Debye ap- 
proximation we obtain 


T ph 


-e(e) = 


47rJV(0)fi D 

m 


Ae 


(18) 


The accuracy of this approximation is known not to be 
high. However it is unlikely to be worse than the use of 
one branch non-dispersive phonon spectrum model that 
we adopted at the beginning. The use of expression (18) 

yields r~} (e) = r ph - e (Q D ) — . In the same approxi- 

^ 1d 3 

mation we obtain from (8) T e _ p / ( (e) = t s ( — — ) . Com- 

\“D/ 


bining the two results we obtain 


Hi — H/} Y T s j Tph—e (I^£)) (19) 

This expression with T p h- e {^D) taken from (18) turns 
out to be close to the one obtained directly from 
T ph—e ( n i) = T e-ph (Hi) and using (4) and (8). Then as- 
suming Hi < Hd and taking for a 2 (H) its limiting value 
a 2 (0) and approximating E(H) by its Debye value, which 
is accurate in this frequency range, we finally obtain 


fir = n D ^®^( l + A) (20) 

This is a better approximation than (19), because the 
only assumption used for its derivation is a 2 (Hi) = 
a 2 (0). For Au we obtain Hi/H£>= 0.036 and 0.039 from 
(19) and (20) respectively. For Mo we obtain from (20) 
Hi/Hii=0.12, which sets a higher threshold than Hi in 
gold. 

Using for an estimate the numbers for Au as 
weaker coupled e-ph system, Ep = 5.51 eV, flp = 


14.2 meV, r s = 3.0 and r s = 130 fs, we obtain 
E*/Hu « 8.2 neglecting the last term in (17) and 
Ei/El ss 4.0[a _1 (l + T s /Tj mp )] 3 / 5 ~ 4.0. The duration 
of Ei — ^ H p stage therefore is tdc <0.8ps. 


B. Modeling the bilayer 

The model parameters determining the phonon loss 
and spectral line-shape are fii/H^, L(E), D } t dc , f3 = 
l P h e (Ho ) / d and r] e ff- The first group are the parameters 
depending on the bulk properties of sensor. The diffu- 
sion coefficient D of down-converting electrons and t dc 
enter in combination Dt dc as a single parameter defin- 
ing the size of excited volume filled with first genera- 
tion of phonons. During the E-[ — > Qp> down-conversion 
stage the diffusion coefficient of a bilayer is modified by 
strong elastic scattering due to phonon emission, so that 
D = D(0)/(l + Tim P /T s ), where -D(O) is the measured low 
temperature value. The X-ray absorption depth L(E) 
and correspondingly distribution function of absorption 
sites P s (zq,E) is modeled in a piece-wise way to incor- 
porate the discontinuity of X-ray absorption in a bilayer. 
This is necessary to account for the factor ~2.5 differ- 
ence in X-ray absorption efficiency in Au and Mo. The 
parameter ij e ff is of special importance. It depends on 
bulk properties such as the acoustic impedances of both 
the sensor film and the substrate. It also depends on the 
properties of the sensor /substrate interface. 


1. Ideal interface 

For ideal atomically clean interfaces without phonon 
scattering, 77 and p and £ c can be calculated using the 
acoustic mismatch model (AMM) of the elastic isotropic 
continuum 35 . This model allows us to evaluate the av- 
erage transmission coefficients for phonons of individual 
polarizations which are incident from the metallic side 
onto the horizontal interface with the substrate. There 
are two shear waves with horizontal and vertical dis- 
placements, SH and SV respectively, and a longitudi- 
nal L-wave. AMM coefficients are useful in defining the 
probable range for phonon transmission of longer wave- 
length phonons depending upon the elastic impedances 
of the metal and substrate and the angles of total in- 
ternal reflection for individual modes. We will use the 
one-mode model of the phonon spectrum, and take the 
sound velocity c of the metal to be equal to the one en- 
tering the Debye model, i.e 3/c 3 = 2/cf + 1/cf , where c t 
and Ci are sound velocities for transverse and longitudi- 
nal phonons respectively. Similarly, for the substrate we 
use c s which is expressed using the same relation through 
respective sound velocities c s t and c s j. Furthermore we 
assume that c is significantly smaller than c s , so that 
the angle of total internal reflection 9 C = arcsin(c/c s ) is 
small. Under these conditions we approximate 77 as for 
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SH waves 35 


and p = sin 2 


Oc 
2 ' 


We stress that the actual 


product pp( 1 + £ c )/2 entering the energy loss and affect- 
ing the spectral lineshape is a fitting parameter to be 
found from fitting the theoretical Ei oss and line-shape to 
experiment, thus such an estimate serves as a consistency 
check. 


2. Rough interface 

The transmission of higher-frequency phonons re- 
quires a special analysis. The first generation of phonons 
peaks at the maximum of phonon density of states. Thus, 
the characteristic energy is ~ f Id- The transmission of 
these phonons may significantly differ from the AMM 
predictions. For a sufficiently rough interface phonon 
transmission is better approximated by the diffuse mis- 
match model (DMM) 36 . In fact real transmission lies 
between the two estimates 37,38 . In the case of roughness 
scattering the expression (7) also can be used, but the 
factor 77p(l+^ c )/2 must be replaced by an effective trans- 
mission coefficient p e ff- A significant difference between 
the two limiting cases arises because of different physical 
mechanisms of phonon transmission. For an ideal inter- 
face the AMM transmission properties are determined 
by the bulk characteristics of the sensor and adjacent 
medii, i.e acoustic impedances. For scattering-assisted 
phonon transfer across the interface the transmission be- 
comes strongly dependent on the correlation scales of the 
interface roughness. Therefore, the transmission charac- 
teristics for first generation phonons and longer wave- 
length athermal phonons of subsequent generations may 
strongly differ. The most important distinction is not 
only in quantitative differences resulting in difference in 
magnitudes of transmission but also that in scattering as- 
sisted phonon transmission of short wavelength phonons 
the interface roughness is being probed locally and thus 
the transmission coefficient may fluctuate from point to 
point along the interface. Thus, Ei oss ^ei^, n D will con- 
tain a stochastic component dependent on lateral coor- 
dinates and introducing an extra broadening to spectral 
line-shape. 


C. “Effective medium” model 

The simplest model for considering phonon trans- 
port in a bilayer is that of the “effective medium”, a 
single layer, with thickness d the same as the thickness 
of Mo/Au bilayer but with phonon characteristics being 
intermediate between the values for Au and Mo. For ex- 
ample, we may let c to vary between 1.4 • 10 5 < c < 
3.7 • 10 5 cm/s, the lower limit being the average sound 
velocity for Au, and the upper - for Mo. The athermal 
phonons entering from Au into Mo layer interact with 
electrons in Mo stronger than in Au so that phonons can 
further down-convert even in a small thickness layer of 


Mo. The down-conversion process in Mo is no different 
to the one described above for Au, and hence we may also 
discuss down-conversion in effective medium in terms of 
effective parameters having intermediate values between 
those in Au and Mo. There exists a slight problem con- 
nected with small number of cases of X-ray absorptions in 
Mo. Phonons of the first generation in Mo above f Id,Au 
in gold (we added extra reference to Au to emphasize that 
this quantity relates to gold) cannot enter and propagate 
in Au and hence cannot be described in the framework of 
the effective medium model. However, these phonons in 
Mo (~ 33 meV) decay very fast and may be disregarded 
on their own. Their decay products, the subsequent gen- 
erations of athermal phonons in Mo, fall below £Id,Au 
and may be considered within the model of an effective 
homogeneous single layer. 


D. “Effective interface” model 


An alternative model can be also suggested. In this 
model we completely ignore X-ray absorptions in the 
Mo layer for the Au/330nm/Mo/50nm bilayer. The ab- 
sorption coefficient in Mo at 6 keV is a factor of 2.4 
smaller than in Au. With thickness of Mo layer being 
~ 0.13d only 6.9% of the X-ray photons are absorbed in 
the Mo. The role of Mo layer in this model is in modi- 
fying the Au/Si transmission characteristics of phonons. 
In this sense the Mo layer plays the role of an “effec- 
tive interface” . The properties of an “effective interface” 
are determined by strong e-ph interaction in Mo. From 
(18) we obtain for Mo layer T. p h- e (^D) = 1-7 ps and 
Tph-eiS^i) — 15 ps. Therefore the mean free paths for 
phonons with f l ~ Qd,Au are short in comparison with 
the thickness of the Mo layer and these phonons entering 
Mo from Au layer are re-absorbed and further experience 
enhanced down-conversion in Mo. A rough estimate for 
the effective transmission coefficient to be used instead 
of rjp( 1 + £)/2 in expression (7) can be taken as 

VAu /M oPAu /Mo ( 1 + £ Mo/Au )/2 . , 

Veff = -i — 1Z ( 21 ) 

i 'esc, Mo/ Si/ 'esc, Mo/ Au 

This expression reflects the balance of the flux incoming 
from Au into Mo with the outgoing phonon fluxes from 
Mo into gold and Mo into the Si substrate. The ratio of 
the latter is taken as the ratio of phonon escape times 
from the Mo layer of thickness dMo via the respective 
interfaces, T esCtMo /Au = 4d M o/r)Mo/AuCMo and similarly 
for the escape into the Si substrate, and is equal to the 
ratio p M o/Au/ VMo/ Si- 


V. MODELING THE SPECTRAL LINE-SHAPE 
IN MO/AU SENSOR FOR MnK a -COMPLEX 

The energy loss from a thin film is given by the ex- 
pression (10). If the source of X-rays is mono-energetic at 
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1-E/fi 1-EEo 

FIG. 7: Two calculated spectral redistributions, given in 
terms of probability versus energy, due to the energy loss 
in a thin film from a mono-energetic photon source. The 
“double-peak” line-shape (arbitrary units) depends upon the 
film material parameters, and two examples are shown. 


Ts PS 

^ph — ei^D^) PS 

c • 10 5 cm/s 

a eV 

Veff 

D cm 2 /s 

tdc ps 

0.13 

80 

1.4 

4.15 

0.035 

500 

0.75 


TABLE I: Model parameters for Au film 


S(E) 




E [eV] 


FIG. 8: (Color online) Simulated normalized line-shape in- 
cluding the contributions from absorption within different 
depth regions of the Au film. 


a photon energy Eq, then the line-shape is determined by 
the coordinate dependent part of Ei oss (zo). As was deter- 

dEi oss (z ) 


mined from expressions (7) and (9) 


dz 


= 0 . 


±d / 2 


This is a consequence of the phonon distribution being 
related to the local density of diffusing electrons which 
are characterized by zero flux at the interfaces with vac- 
uum or dielectric. For a — > 0 we have S(z,Eq) —>■ 
6{E-E 0 + Ei oss (z)) = 6{z-z*)/\E' loss (z*)\, where z*(E) 
is the solution of E — E 0 + Ei oss (z) = 0. Due to the pres- 
ence of derivative in the denominator there appears an 
integrable singularity at E = E 0 + Ei oss (±d/ 2) resulting 
in a characteristic “double-peak” or shoulder in the line- 
shape as depicted in Fig. 7. The left and right panels 
correspond to different material parameters. 

We will model our experimental results using the 
Holzer spectrum 39 for the MnK a -complex. The “double- 
peak” structure will be present for all components of the 
Holzer spectrum. Therefore, the resulting line-shape is 
more complicated and contains more features. These fea- 
tures however may be significantly smeared if the baseline 
resolution described by the variance er is not sufficiently 
high. Below we describe in detail simulation of the line- 
shape using the “effective interface” model to describe 
Mo/Au bilayer. Table I shows all parameters used for 
the line-shape simulation. 

The choice of parameters has been discussed in the 
previous section. For rj e ff in Table I we show the AMM 
estimate for Au/Mo interface reduced by the fitting fac- 
tor 2 (see below) to account for the role of the denomi- 
nator in (21). Fig. 8 shows the results of the simulation 
together with a breakdown of the overall spectrum into 
contributions from different depths of the 330 nm Au 
film, i.e. from four different notional 82.5 nm thick sub- 
layers. The outmost left line is for the sub-layer which 



FIG. 9: (Color online) Comparison of the simulated line and 
experiment 


is neighboring the Mo layer and shows the highest loss. 
This is followed by successive sub-layers at larger dis- 
tances from the escape interface and therefore smaller 
respective losses. Fig. 9 shows the simulated line super- 
imposed upon the experimental result. 

When we compare our simulated curve with the ex- 
perimental data it is apparent that both have approx- 
imately the same FWHM. Thus the dependence of en- 
ergy loss on the coordinate of absorption site results in 
line broadening which is roughly in agreement with the 
experiment. The simulated line-shape shows a signifi- 
cantly different behavior on both sides of the peak and 
also a well defined shoulder. The origin of the shoulder 
was discussed above and is clear from Fig. 8. The ex- 
perimental line exhibits a longer tail on the lower energy 
side and also a much smoother drop on the left side of the 
peak. The theoretical curve also exhibits a rapid decrease 
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on the right-hand side of the peak, with this drop occur- 
ring over 10-15 eV, in contrast to the measured spectrum 
which has a much slower decrease in number of counts 
with energy. 

The lower energy tail is easier to interpret within the 
model of energy loss from Au/Mo him. Absorptions in 
Mo typically results in a bigger loss from the Mo him. 
This is not only due to the closeness to the escape inter- 
face, but also due to the better acoustic matching to the 
Si substrate. The additional broadening to lower energies 
in Fig. 9 can also be explained because the distribution 
function P s {zq, E) for a bilayer has not yet properly been 
taken into account. It can be written in the form 


Y Au exp(— T^d) 


Ps ^ z ° ,E ^ l - exp(-T AudAu - T ModMo) 
[eKp(-T Au z 0 )Q(d/2 - z 0 )Q{z 0 + d/2 - d Mo )+ 


exp(— r Mo 2 O ) 0 (z o + d/2)Q(d Mo 

-L Mo 


d/2 - zq) 


(22) 


accounting for the weaker absorption in Mo relative to 
Au. This will reduce the weight of the smallest (left) 
peak in Fig. 8 ( due to absorptions close to the escape 
interface), down by a factor of 2.4. 

It is not possible to explain the higher energy side 
of the peak profile within any of the two models. With 
the experimental baseline resolution a = 4.15 eV and the 
intrinsic distribution of Mn K Q X-rays (as described by 
Holzer) extending over just 10 eV, there are no physical 
reasons within our model of phonon losses to have this 
profile less steep than 10-15 eV. 

In order to understand the leading edge discrepancy 
we have examined the phonon escape from the sensor in 
more detail. As we have seen, the spectral line-shape is 
sensitive to both the bulk phonon transport character- 
istics and the transmission across interfaces. Simplifica- 
tions in modeling the bulk phonon transport cannot be 
responsible for the failure to correctly reproduce all the 
major features of the spectral line-shape. More realistic 
models of phonon spectrum and evolution of phonon dis- 
tributions will only slightly modify the statement that 
the first generation of down-conversion phonons has a 
dominant effect on spectral line-shape, while subsequent 
generations due to their fast spatial spread only con- 
tribute to the height of the peak and its horizontal shift 
(total energy loss). Our derivation of the distribution 
function for the first generation of phonons can be gen- 
eralized to a spectrum containing 3 acoustic branches. 
A generic feature of this distribution function is that it 
is formed as a result of spontaneous emission of phonons 
by down-converting electrons during the — >• YId stage, 

when electrons are capable of emitting all phonons. It is 
also important that the lifetimes of all phonons emit- 
ted are much larger than the duration of this stage. As 
a result, the first phonon generation distribution func- 
tion is a snap-shot of distribution of phonons at a time 
t = tdc- The details of this distribution are determined 
by the densities of states at each energy together with 


the Eliashberg coupling functions for individual modes 
to determine the relative magnitudes. The more compli- 
cated three branch bulk model is capable of changing the 
structure on the left-hand side of the measured X-ray 
spectrum by splitting the shoulder into individual fea- 
tures due to loss of transverse and longitudinal phonons. 
However, it is unlikely that this structure will be smooth 
and it definitely cannot introduce a smoother variation 
for the higher energy side of the X-ray spectrum. On 
the left, the spectrum begins from X-ray absorption from 
E = E 0 — Ei oss (—d/ 2), while on the right it terminates 
at E = E 0 - Ei oss (d/ 2). 

The origin of the discrepancy between the model 
and data is most likely to be connected with the trans- 
mission characteristics of phonons of the first generation 
whose escape determines the spectral line-shape. So far 
we have shown the results of simulations in Figs. 8 and 
9 assuming coherent transmission of phonons across the 
ideal interface. Phonons of the first generation in Au 
have energies close to Hd, corresponding to maxima in 
the phonon density of states in Au of 6-8 meV for trans- 
verse and 12-14 meV for longitudinal phonons. The wave- 
lengths of these phonons are on the nanometer scale so 
that any roughness of the interface on the nanometer 
scale will cause strong phonon scattering. The transmis- 
sion of the first generation phonons < which are 

incident onto the Au/Mo interface from gold, is more 
strongly affected by the roughness than the transmission 
of the same frequency (but larger wavelength) phonons 
which are incident onto the same interface from molyb- 
denum side. Moreover in the “effective interface” model 
of Mo/Au bilayer phonons inside Mo layer experience 
more effective down-conversion and therefore at any mo- 
ment of time will have lower mean frequencies. Thus, the 
phonon transmission across the Au/Mo interface (from 
Au to Mo) may be determined by its roughness while at 
the same time transmission of phonons across the Mo/Au 
interface, i.e. moving in the opposite direction (from Mo 
to Au) can be described by the AMM model. If this is the 
correct picture then the expression (21) can be re-written 
as 


Vs,Au/ Mo /V Au/Mo 
Veff = V Au/Mo i . 7 

-L r VMo/Au/'Hesc,Mo/Si 


vriAu/Mo (23) 


where T] s , Au/Mo is the scattering assisted transmission co- 
efficient for high frequency phonons from Au into Mo, 
while ij A u/Mo i V Mo/Au and VMo/Si are AMM transmission 
coefficients and v is the fitting parameter. The scatter- 
ing assisted effective transmission coefficient, rj s ^ Au/Mo , 
differs from its AMM counterpart. Correspondingly, v 
is likely to be less than 1 because of the large value of 
the ratio r/Mo/Au/Vesc, Mo/Si being balanced by the fac- 


?/s, Au/Mo 

tor 


> 1. In our simulation of the line-shape in 


V Au/Mo 

Figs. 8 and 9 we chose v = 1/2. The coupling constant 
entering the scattering cross-sections may be large due 
to strong variation of material densities at the interface. 
An important aspect of scattering assisted phonon 
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FIG. 10: (Color online) Comparison of simulation with exper- 
iment. Averaging the pulse in Fig. 9 over lateral coordinates 
with a v = 0.065 fj e ff 

transfer is the statistical variation in the probability of 
scattering and the dependence of the scattering cross- 
section of short wavelength phonons on the local char- 
acteristics of the roughness of the interface. This intro- 
duces statistical fluctuations in the interface scattering 
and hence in the transmission coefficient in longitudi- 
nal directions. Thus, me may expect the effective trans- 
mission coefficient to be a random function of the lat- 
eral coordinate r, r/ e //( r). We have implemented the 
averaging over the plane of the interface assuming a 
Gaussian distribution of r] e f / with the variance a v about 

mean value fj eff , P v = exp [-(%// - ?7 e //) 2 /2cr^] . 

v r\ 

Fig. 10 shows the result of this simulation. Overall the 
agreement between the model and data looks much bet- 
ter. The slopes of the leading and trailing edges of the 
distribution are reproduced well with some tailing in the 
data still remaining on both sides. 

Some of the remaining tailing can be explained by 
the model not having taken into account that 6.9% of 
X-ray absorptions generate phonons in the molybdenum 
part of the Mo/Au bilayer. The enhanced loss from Mo 
can result in a small separate peak on the left side of the 
spectrum. We could implement this by attributing some 
extra energy loss for all absorption events in the layer 
closest to Si interface. In Fig. 8 these events are from 
the notional sub-layer on the left. By including some 
small extra energy loss for absorption in this sub-layer 
the corresponding spectrum is shifted further to the left, 
potentially giving a better fit to the data. 

The broadening seen in the leading edge of the data 
in the region 3 can be due to the fact that effective trans- 
mission has both coherent and scattering mediated re- 
gions. The presence of disorder-free patches of the in- 
terface, where this region is a small fraction of the total 
interface area, will mean that 7y e // is locally suppressed, 
and will result in smaller energy loss. Some of the re- 
maining tailing in region 1 is difficult to explain within 


the phonon loss model and probably has a different ori- 
gin. 

The extended tail in the region 1 can be interpreted 
as being due to universal mechanism related to loss of en- 
ergy during the E 0 — >• E\ cascade. A certain percentage 
of the photon energy can leak via the loss of energetic 
electrons if their energy is larger than the work function 
from a metal by the time they hit vacuum surface or di- 
electric interface. Our modeling of this mechanism of loss 
will be published in separate paper. 

Finally, simulation based on “effective medium” 
model yields very similar results using same parame- 
ters as in Table I except taking for mean sound velocity 
c = 2.0 • 10 5 cm/s. This value is close to weighted value 
for Au/Mo bilayer accounting fractional volumes of Au 
and Mo. 


A. Energy loss and e-e interaction 

So far we have not included in the modeling of 
the spectral line-shape the coordinate independent term 
Eioss , Incorporation of this term results in shift- 
ing the pulse towards lower channels without changing 
its shape. The only extra fitting parameter entering 
this term as discussed in Section III is the threshold 
G-i ■ Including small energy loss, E iosS! thermal , due to 
the additional thermal conductance of heat out of the 
microcalorimeter due to the MPT film being at a higher 
temperature for a short duration as discussed in the ex- 
perimental section we obtain for the total energy loss 
Eioss — Pioss, thermal T E[ oss i -f- Ei oss q d —^q 1 . For the 
set of parameters in Table I the averaged loss due to 
phonons of the first generation is Ei OSSt i = 0.7%. Taking 
Eioss, thermal = 0.4% as estimated above we obtain that 
in order to fit experimental result of Ei oss ~ 3.8% the 
remaining loss term must account for 2.7% of the energy 
loss. Thus, as expected, this is the dominant loss term. 
Fitting this number yields the threshold fli by factor 1.24 
higher than given by expression (20) derived for the sit- 
uation when T ee (n i) T p h~e( Hi). We therefore arrive 
at the conclusion that the electron-electron interaction in 
our Au film is relatively weak and only slightly modifies 
the lower threshold of phonon stage of down-conversion 
cascade. 

This result further proves the consistency of our 
model. Indeed, at higher threshold E\ the estimated 
electron-electron relaxation time was T ee (Ei) — 65 fs. 
We have shown in Section III that electron-electron in- 
teraction rate scales from this energy down practically to 
the fli threshold according to e 3 / 2 law. Accounting the 
correction of 1.24 determined from fitting energy loss to 
experiment we estimate ratio Ei/fli ~ 675. Thus, the e- 
e interaction time scaled down to energy fli is estimated 
to be T ee (f2i) ~ 1.1 ns being close to our initial estimate 
Tph-eiSh) — 1-5 ns. The latter estimate was obtained by 
taking r ph _ e [p, D ) = 80 ps from Table I and scaling its 
rate to the threshold fli according to linear law (18). 
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T s PS 

/ T~ph—e{^‘D^) PS 

c • 10 5 cm/s 

V eff 

D cm 2 /s 

tdc ps 

0.0125 

13.3 

3.18 

0.043 

35 

0.40 


TABLE II: Model parameters for W film 


One of the important practical implications of this 
result is that e-e interaction in 330 nm thick Au film has 
been found to be relatively weak allowing therefore for an 
extended long lasting phonon stage of down-conversion 
cascade. It is this property of high quality Au him of 
Mo/Au bilayer that made possible observation and iden- 
tification of athermal phonon loss. 

VI. DISCUSSION 

Strong energy leak from a thin him was experimen- 
tally observed hrst by Cabrera 40 in a thin him tungsten 
TES designed for detection of single optical photons. The 
measured fraction of energy detected with a W TES was 
surprisingly low ~ 42%. This was attributed to the loss 
of athermal phonons. Such a big energy loss and re- 
duction in energy efficiency is well known for detectors 
operating in non-equilibrium regime, for example, STJs 
and MKIDs. These sensors operate at low temperatures 
T<T C . Generation of non pair-breaking phonons during 
the down-conversion cascade in a superconductor with fi- 
nite gap A results in approximately 42% energy loss due 
to the de-coupling of non-pair-breaking phonons from the 
condensate and the very weak interaction with quasipar- 
ticles. In a TES which is operated close to the critical 
temperature T — T c <C T c the order parameter A — > 0, 
and all the non-equilibrium phonons generated can be re- 
absorbed by quasiparticles. This situation is very similar 
to the electron-phonon interaction in normal metals. 

We are able to explain this high loss of energy using 
the same calculations used to explain the 3.4% energy loss 
from the escape of athermal phonons in our Mo/Au MPT 
bilayer. On this basis, we suggest, for the first time, that 
the large energy loss from a W TES can be explained 
by athermal phonon loss when taking into account all 
possible modes of athermal phonon loss. In an earlier 
analysis of spectral degradation of optical TESs due to 
phonon down-conversion noise 4 the energy efficiency in 
W TES described by Cabrera et al was estimated to be 
7%. However, this calculation only took into account 
the first two phonon generations. It is now clear that a 
dominant contribution to energy loss can be due to lower 
energy phonons which are emitted during the extended 
fi/) — >■ fli stage. 

Using the calculated Eliashberg function for 
tungsten 41 and parameters from Table II after averag- 
ing over the coordinates of absorption cites we obtain 
from (5), (7), (9) and (10) Ei oss ~46.9%. The break- 
down of the observed losses due to different mechanisms 
for this experiment requires a special analysis. However, 
the overall number that we have calculated clearly domi- 


Ts ps 

Tph—eiS^'D') PS 

c ■ 10 5 cm/s 

creV 

D cm 2 /s 

tdc ps 

0.026 

1.2 

2.33 

3.8 

500 

0.25 


TABLE III: Model parameters for Ta film 


nates the observed 58% loss. It is also consistent with the 
magnitude of the experimentally observed improvement 
of energy efficiency from ~ 42% to ~ 80% that was seen 
in TESs that were fabricated on micromachined SUN4 
membranes 1 , which reduced the athermal phonon loss. 
The scattering rates derived from the parameters shown 
in Table II are consistent with the estimates of electron- 
phonon, phonon-electron and electron-electron scattering 
rates in tungsten 42 . Both sets of data give close values 
for fli and r e _ p / l (fl 1 ) <C T ee (Oi) justifying our estimate 
energy loss. 

It is remarkable that a relatively simple model can 
correctly describe the loss of athermal phonons in the 
process of energy deposition from such different thin 
metal films, even when the resulting fractional energy 
loss efficiency varies by over an order of magnitude. 

In an experiment using an STJ consisting of 
Ta(110nm)/Al(5nm) /A10x/Al(5nm)/Ta(95nm) grown 
on a sapphire substrate by Verhoeve et al. 43 , the best 
spectral resolution achieved was A E = 15.7 eV [FWHM] 
at 6 keV. Based upon the responsivity of the STJ, and 
the total noise from all noise contributions, the best 
achievable energy resolution was 8.9 eV [FWHM]. 

In Fig. 11 we show the results of spectral shape mod- 
eling for the geometry of this STJ. The material parame- 
ters for Ta used for this modeling are shown in Table III. 
Just as we did for the previous experiments, the intrinsic 
line-shape of MnK a X-rays from Holzer was used 39 , and 
a detector resolution of 8.9 eV [FWHM] was assumed. 

In order to evaluate t s we used the calculated Eliash- 
berg function and its first reciprocal moment A for Ta 34 . 
The results from this modeling is shown in Fig. 11, as 
well as the original spectrum. It is apparent that the 
co-ordinate dependent energy loss via escape of athermal 
phonons can adequately explain most of the measured 
spectrum. 

In another similar experiment using a Ta-based X- 
ray absorber on a solid substrate, utilizing a distributed 
STJ detector read-out, Li et al. 44 reported a spectral res- 
olution of 12 eV [FWHM] for 6 keV X-rays, close to the 
theoretically achievable limit for this particular device. 
According to our models, the expected broadening that 
occurs due to the ^-dependent athermal phonon energy 
loss for this 600 nm thick Ta absorber is ~ 6 eV [FWHM]'. 
Therefore this level of spectral broadening would not sig- 
nificantly impact the performance of this sensor. 

VII. CONCLUSIONS 

We have shown experimentally that following a pho- 
ton absorption in a thin metal film grown on a solid 
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Channel 


FIG. 11: This graph shows the spectrum of MnK a observed 
by Verhoeve et al . 43 The solid line superimposed on this spec- 
trum is the result of our simulated line-shape which takes into 
account the spatial broadening of MnK„ complex due to loss 
of athermal phonons as well as intrinsic and electronic noise 
sources. 


experimental line-shape allows us to evaluate the degree 
of scattering-assisted phonon transmission due to sub- 
nanometer-scale roughness of the interface. 

Our relatively simple model explains a wide variety 
of measured energy losses and effects on measured spec- 
tral resolution. It is important for the design of future 
detectors with small or large regions of solid substrate 
contact. 
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substrate, a substantial fraction of the energy deposited 
may be lost prior to establishing thermal equilibration 
due to the escape of athermal phonons to the substrate 
below. This form of energy loss can significantly affect 
the performance of high-resolution spectrometers on solid 
substrates. When the whole energy of the X-ray was re- 
tained in the spectrometer following its absorption in the 
absorber the measured spectral resolution of MnK Q com- 
plex was 4.3 eV [FWHM]. However, for absorption events 
directly in the sensor film with measured energy loss 3.4% 
the energy resolution deteriorated to « 40 eV [FWHM], 
The main source of degradation in soft X-ray region is 
now understood to be due to the “z” -dependence energy 
of athermal phonon loss, while the effect of fluctuations 
in energy loss is small. 

We have developed a theoretical model for simu- 
lation of the spectral line-shape and have shown that 
the latter is sensitive to details of energy deposition, 
energy exchange between electron and phonon systems 
and phonon transmission to the substrate across the es- 
cape interface. Apart from trivial geometrical factors the 
amount of energy loss reflects: the relative strength of 
electron-electron versus electron-phonon interaction over 
a broad spectral range, ~ 1 eV, about the Fermi en- 
ergy and the magnitude of the phonon transmission co- 
efficient. The spectral line-shape is mostly determined 
by the high energy, short (atomic scale) wavelength first 
generation of phonons which are released during the ini- 
tial stages of the energy down-conversion process. The 
loss of lower-energy athermal phonons of subsequent gen- 
erations results only in a parallel shift of the spectral line 
without affecting its shape. Fitting the theoretical to the 
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